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Abstract 

An analytical study of discrete and resonant impurity quasi-Coulomb 
states in a narrow gaped armchair graphene nanoribbon (GNR) is per- 
formed. We employ the adiabatic approximation assuming that the mo- 
tions parallel ("slow") and perpendicular ("fast") to the boundaries of 
the ribbon are separated adiabatically. The energy spectrum comprises a 
sequence of series of quasi- Rydberg levels relevant to the "slow" motion 
adjacent from the low energies to the size-quantized levels associated with 
the "fast" motion. Only the series attributed to the ground size-quantized 
sub-band is really discrete, while others corresponding to the excited sub- 
bands consist of quasi-discrete (Fano resonant) levels of non-zero energetic 
widths, caused by the coupling with the states of the continuous spectrum 
branching from the low lying sub-bands. In the two- and three-subband 
approximation the spectrum of the complex energies of the impurity elec- 
tron is derived in an explicit form. Narrowing the GNR leads to an in- 
crease of the binding energy and the resonant width both induced by the 
finite width of the ribbon. Displacing the impurity centre from the mid- 
point of the GNR causes the binding energy to decrease while the resonant 
width of the first excited Rydberg series increases. As for the second ex- 
cited series their widths become narrower with the shift of the impurity. 
A successful comparison of our analytical results with those obtained by 
other theoretical and experimental methods is presented. Estimates of 
the binding energies and the resonant widths taken for the parameters of 
typical GNRs show that not only the strictly discrete but also the some 
resonant states are quite stable and could be studied experimentally in 
doped GNRs. 
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1 Introduction 



The electron properties of two-dimensional (2D) graphene, a single-layer carbon 
sheet, has attracted much attention by both theoreticians and experimentalists 
(see [T] and references therein). Along with this related structures, namely 
graphene nanoribbons, are also under intensive investigation 0. One of the 
reason for this is that the long electron mean free path in graphene up to 1 /im 
opens a field of carbon-based nanoelectronics, where GNRs are used as inter- 
connects in nanodevices. The unique electron mobility in graphene structures 
is caused by the strong bonding between the carbon atoms in the honeycomb 
lattice of graphene. This in turn prevents the replacing of the carbon atoms by 
alien ones. Nevertheless, graphene is not immune to extrinsic disorder and its 
transport properties [3] are very sensitive to impurities and defects [4 . 

The theoretical problem of an impurity in 2D graphene was considered orig- 
inally in O [6l [21 m E] . In the vicinity of the Dirac points in k space, which 
are peaks of the double cones of the Fermi surface, the low-energy electronic 
excitations in gapless graphene are described by the equation of the effective 
mass approximation, which is formally identical to the 2D Dirac equation for 
a massless neutrino, having the Fermi speed vp — lO^m/c. In the presence of 
an attractive impurity centre of charge Z screened by a medium of the effective 
dielectric constant e^Q the electron states are drastically different for the sub- 
critical Jc < J and supercritical Jc> J regions of the strength of the Coulomb 
interaction, where J = \j\ = 1,2, .. . and Jc = Ze^ / AireQe^^F are the 2D 
momentum of the impurity electron and that, having the speed VF^ respectively. 
Clearly, the super- and subcritical regimes can be reached for the dimensionless 
Coulomb potential strength q — Ze^ / ^-neQe^QTiVF for q > |j| and q < \j\, re- 
spectively. The difference of the subcritical and supercritical electron states is 
caused by their different behaviour in the vicinity of the impurity centre f — >■ 0. 
The subcritical regime admits regular solutions to the Dirac equations, while 
the wave functions corresponding to the supercritical case oscillate and do not 
have any definite limit. The physical reason for this is that at the subcritical 
strengths g < \j\ the centrifugal potential barrier prevents the electron "fall to 
the centre" [10], while the supercritical strengths g > \j\ provide the collapse. 
Clearly, the continuum approach based on the Dirac formalism becomes inappli- 
cable. The lattice-scale physics dominates that in turn requires a regularization 
procedure, namely the cutoff of the Coulomb potential at short distances tq ~ a, 
where a ~ 1.42 A is the C-C distance in graphene. The physics of the supercrit- 
ical impurity electron in graphene ,6j closely resembles that of the relativistic 
electron in an atom having the nuclear charge Z > 137 [TTl [T^l [131 [S] . Since, as 
it follows from below, only the supercritical regime is relevant to the impurity 
state in CNR, we focus on this case. 

Numerical and analytical approaches developed on the tight-binding model 
of the graphene lattice and of the Dirac equation, subject to the regularization 
procedure, respectively, undertaken originally by Pereira [7] et.al. have revealed 
the infinite number of the quasi-bound states, having the finite width, arising, 
as it was shown quasi-classically [9], from the collapsed states. If the require- 
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ment of the regularity of the wave functions in the vicinity of the source of the 
electron attraction is to be replaced by the less rigorous condition of its square 
integrability an infinite number of the strictly discrete energy levels were found 
to occur. The Coulomb potential cutoff in the gapped |T3] and no cutoff in the 
gapless graphene [IS] induce the energy series bounded and unbounded from 
below, respectively. 

In the GNR, which in principle can be treated as a quasi-lD structure, we can 
expect completely different results. The strictly discrete bound states regular at 
the impurity centre (r = 0) are realized without the regularization procedure, in 
particular without the cutoff of the Coulomb potential, preventing the collapse. 
Consequential concerns lie in the well known fact that the reduction of the 
dimension of the structure increases the stability of the impurity electron. In 
units of the impurity Rydberg constant Ry the binding energy Ef, of the impurity 
electron in 3D bulk material is i^f, = 1, in the narrow 2D quantum well Eb — 4: 
[T7] , and in the thin quantum wire of radius R much less than the impurity 
Bohr radius oq, -Eb ~ ln^(i?/ao) [IH]. Besides, an extremely weak 3D atomic 
potential not providing bound electron states, transforms in the presence of a 
magnetic field into a quasi- ID system binding the electron [19]. It is relevant to 
note that these atomic states arise under as weak as one likes magnetic fields i.e. 
the as large as one likes magnetic lengths playing the same role as the width of 
the GNR. Note that the confinement attributed to the semiconductor thin films 
[5D] replaces the 3D Coulomb potential r~^) by the effective 2D potential of 
the weaker singularity of the logarithmic character. It seems that in the quasi- 
lD GNR the effect of the attenuation of the potential singularity preventing the 
fall to the centre exceeds that of the vanishing of the 2D centrifugal potential 
barrier promoting the collapse. 

Clearly a study of the impurity electron state in graphene structures is im- 
portant on account of two reasons. First, these structures provide a realization 
in solid state physics of remarkable effects of quantum electrodynamics caused 
by a large "fine structure constant" /Tivp ~ 2.5 [5] HI HH [22]- Second, we 
expect a strong impact of impurities on the electronic systems not only for 2D 
graphene layers possessing an outstanding high electron mobility [6] but in par- 
ticular for impurity GNRs whose properties are not widely addressed in the 
literature yet. 

Brey and Fertig [J^ have shown that the energy spectrum of the electron 
in an armchair GNR bounded in x-direction is the sequence of the subbands 
formed by the branches of the continuous energies of the longitudinal unbounded 
j/-motion emanating from the size-quantized energy levels , {N is the discrete 
lable), reflecting the ribbons a;-confinement. The equation for the components 
UA.B{y) of the Dirac spinor relevant to the A and B sublattices of the graphene 
for the electron positioned far away from the impurity centre has the form 

h vj, 

showing that the armchair GNR manifests itself as gapped structure entailing 
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Figure 1: A schematic form of the potentials V/vjv(2/) provided in eqs. ([2]), 
(fTT|) . at xo = and quasi-discrete n ([25]) and continuous fc^r spectra 
adjacent to the ground = (discrete states) and first = 1 and second 
A^ — 1 size-quantized levels £jv in the GNR of width d. 



the bound and unbound impurity states for the energies E"^ < e% and E"^ > e%, 
respectively. 

Of special interest is the narrow GNR of width d for which d << rp, where 
rQ{d) ~ h^VplE^ — e^)~^ is the radius of the bound electron state, being in- 
duced by the ribbon confinement d. Such a GNR provides the expected elec- 
tron binding energy Ei, eNf{q), where f{q) is a some function vanishing at 
g = 0, which is of interest and attractive because of two aspects. On the one 
hand the ribbon provides a considerable impurity binding energy which could 
be measured experimentally and on the other hand the impurity potential can 
be treated perturbatively and an analytical approach to the problem becomes 
feasible. 

A comment concerning the from of the energy spectrum is in order. In the 
zeroth approximation of isolated size-quantized A^-subbands i.e. in the single- 
subband approximation the slow longitudinal motion parallel to the bound- 
aries is governed by the 2D Coulomb potential averaged with respect to the 
"fast" transverse A^-states. The energy spectrum consists of sequence of se- 
ries of quasi-Coulomb discrete A^n-levels and continuous sub-bands positioned 
below and above, respectively relatively to the size-quantized energy levels ejv 
(see Fig. 1). Only the series of the impurity energy levels Eon < Eq adjacent 
to the ground size-quantized energy level Eq is strictly discrete. The A^n-series 
adjacent to the excited levels A^ > come into resonance with the states of 
the continuous spectra of lower subbands and in fact in the next multi-subband 
approximation turn into quasi-discrete resonant states (Fano resonances) [24] . 
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The corresponding resonant widths Tj^n determine the aiito-ionization rate and 
hfe-time T^n = Ti/^Nn of the resonant impurity states being of relevance to an 
experimental study. Also to our knowledge an analytical approach to the prob- 
lem of impurities in GNR providing the explicit dependencies of both discrete 
and especially quasi-discrete electron states on the width of the GNR (d) and 
the position of the impurity centre within the GNR are not comprehensively 
available in the literatiire. 

In order to fill the above mentioned gap we perform an analytical study of 
the strictly discrete and resonant impurity states in a narrow armchair GNR. 
The impurity centre is positioned anywhere within the ribbon bound by the im- 
penetrable boundaries. The width of the GNR is assumed to be much less than 
the radius of the impurity state. The complete 2D envelope wave function satis- 
fying the masslcss Dirac equation is expanded with respect to the basis formed 
by the ID size-quantized subband wave functions describing the fast transverse 
motion bound by the boundaries of the GNR. The generated set of equations 
for the ID quasi-Coulomb wave functions relevant to the longitudinal slow mo- 
tion is solved in the single-, two- and three-subband approximations, in which 
the ground, first and second excited subbands are involved. The mathematical 
method is based on the matching of the Coulomb wave functions with those 
obtained by an iteration procedure at any point within the intermediate region 
bound by the ribbon width and the radius of the created Coulomb state. Both 
the real and imaginary parts of the complex energy levels are calculated in a 
single procedure. The dependencies of the binding energy and resonant energy 
shift and width on the width of the GNR and the position of the impurity centre 
are obtained in an explicit form. Numerical estimates show that for a narrow 
GNR the binding energy and the resonant width are quite reasonable, to render 
the impurity electron states in GNR experimentally observable. Our analytical 
results are in line with those calculated numerically and revealed in an experi- 
ment. We remark that our aim is to elucidate the physics of the impurity states 
in GNR by deriving closed form analytical expressions for their properties. We 
do not intend to compete with the results of computational studies. 

The paper is organized as follows. In Section 2 the analytical approach based 
on the multi-subband approximation is described. The real quasi-Coulomb func- 
tions of the discrete and continuous spectriim and the real energy levels deter- 
mining the binding energies are calculated in the single-subband approximation 
in Section 3. The complex energies including the resonant shift and width as- 
sociated to the first and second excited subbands are found in the double and 
three-subband approximation, respectively, in Section 4. In Section 5 we discuss 
the obtained results. Section 6 contains the conclusions. 

According to the above an analytical description of the stable and metastable 
impurity electron states in the narrow armchair GNR is of significant interest. 
It elucidates the underlying basic physics of the carbon-based nanodevices, in 
which the highly mobile electrons remain unbound in the 2D graphene mono- 
layers while in their interconnects, namely in quasi-lD armchair GNRs these 
electrons are trapped by impurity centres. The latter could modify the overall 
transport properties. 
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2 General approach 



We consider a ribbon of widtli d located in the x — y plane and bounded by the 
lines X = ±d/2. The impurity centre of charge Z is displaced from the mid-point 
of the ribbon a; = by the distance —d/2 < xq < d/2. The equation describing 
the impurity electron at a position p = {x,y) possesses the form of a Dirac 
equation 



H(/c, p)^ip) = E^{p); k = -iV; 
where the Hamiltonian H is given by 



H =p 



Ho(A:,) + f^i{ky) + V{p)l- p = hv;v^ 10^ m/c 



with 



(1) 



composed by the Hamiltonians relevant to the inequivalent Dirac points {—ak) 
and (— (T*fc) iff = (<7x,(Ty) are the Pauli matrixes) presented originally in 

Ref . ^31 . The matrix I in ([1]) is the unit matrix and 

V{p) = -^=J===; p^^^- (2) 

is the 2D Coulomb impurity potential, Egg is the effective dielectric constant 
related to the static dielectric constant e of the substrate by [BJ [25] 

The envelope wave four-vector 'J consists of two vectors a.b describing the 
motion of the electron in sublattices A and B of graphene 

^(p) = ^p^(p) + #b(p) ; (3) 

each determined by the wave functions "0^^ ^ 



J Ui'. 

The total A, B state implies the multiplication of the V'a^b "* functions with 
the factors exp{ii?'^"'"'~^/9}, respectively. The boundary conditions for the arm- 
chair ribbon require the total wave function to vanish at both edges x = ±d/2 
for both A, B superlattices [T] 
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(4) 



where — {±K, 0) ; K — 47r/3ao, oq — 2.46 A is the graphene superlattice 

constant. 

The basis wave vectors $ jv (2;) and the energies e n describing the transverse 
size-quantized x-states are derived from equation 



Ho(fc^)'I>Ar(a:) =eAr$(a:) 



to obtain 



where 



with 



<^na{x) + <^NBix) ; $jvA(a; 



(+) 
Vna 





(pmix) 



(+) (-)* (+) 



x-{N-d) 



(— )* 

-'/'TVS = VJVO 



exp < 1 



Kd 



En 



= 0,±1,±2, 



TT 
TT 



(5) 



(6) 



(7) 



We consider transverse states with positive energies En > in the armchair 
ribbon of width d providing the gaped (insulator) structure 7^ [1] . It foUows 
from eq. ([7]) that the energy levels e^v as a function of width d are the oscillations 
describing by parameter a{d) imposed on the decreasing curve ~ d^^. Below 
we ignore these oscillations keeping a — const. 

The boundary conditions (jl]) after substitution of iPa^'b ^ V^na b ' respec- 
tively, are satisfied. The wave vectors ^na,b form orthonormal subsets, for 
which 



(B) 

{^N'A{B)\^NA{B}) 



£N^NB(A) ; {^N'A{B)\^NB{A)) = ; 

Sn'n ; {^n'I^n) = 5n'n- (8) 



The boundary conditions (U)) imposed on the wave vector ^'(p) ([3]) force us 
to expand the wave vectors ^'a,b in series 
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with respect to the basis functions lys^^ si^) taking for the coefficients 



'■NA,B 



(9) 

(y) = 



'■NA.B 



(y) = UNA.Biy)- Substituting the wave vector 5* (|3]) with the wave vectors 

"^A.B and the wave functions "0^3 ■* {p) into eq. ([T]) and subsequently using the 
properties we obtain by the standard method the set of equations for the 
wave functions 



1 / 

—7={UNB + una)] 

v2 



(2) 



1 



—i={uNB - una) ; 
v2 



dy 
dy 



^{E + SN- VNNiy)) v'j^\y) + i Eiv'^iv VN'N{y)v);! (y) 

^AE-SN~VNNiy))v^N\y)~^j: 



VN'N{y)v''N}iy) 



; 
; 



(10) 



VN'Niy) 



V{p) cos 



(N- N')tt 



where the potential V{p) is given by eq. At |?;| ^ d 



dx 



(11) 



VN'Niy) 



\y\ 



Sn' 



N 



o 



3\N'-N\{2s+1) 



; 5 = 0,1,2,.. 



(12) 



y of free electrons in the armchair nanoribbon 



As expected in the limiting case d -> cqs. pUj) decompose into the sets 
describing the ID Coulomb states, while in the absence of the impurity centre 
(Vn'n — 0) we arrive at the wave functions una,b ~ exp(±ifcj,?/) and the 
energies E%{ky) = e% + p^k. 

Below we solve the set (jlOp in the adiabatic approximation implying the 
longitudinal y-motion governed by the quasi-Coulomb potentials VN'Niy) to be 
much slower than the transverse x-motion affected by the boundaries of the 
narrow ribbon. The Coulomb potential ([2]) is assumed to be small compared to 
the ribbon confinement. In the case of ct < 0.5 the lowest three sub-bands are 
specified by indices TV = 0, 1, —1. The set (fTO|) corresponding to these subbands 
becomes 



V 



(2)' 
^0 

(1)' 



+ l{E-eo- Vao) 



(1) 



liE + e,-Vi)v^'^ + ^ 



V 



(2)' 



P 



(1) 



P 

^l.iE~s,-Vi)v[''~l 

p 



(1) 



Vo-ivj^^ 



= 
= 
= 
= 

(2) 



Vi-iv\ 



= 

= 



} (13) 
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3 Single-subband approximation 



At the first stage we neglect the coupling between the states corresponding to 
the subbands with different N. The reason for this is that in the narrow ribbon 
of small width d the diagonal potentials Vnn dominate the off-diagonal terms 
Vn'n {N' ^ N) almost everywhere but for a small region \y\ < d (see eq. (HH)). 
In this case Vn'n — Vnn^n'N and the set ([TU)) decomposes into independent 
subsets each specified by an index N. The ID impurity states are then governed 
by the potential 




di^2 = d±2xo ; < xq < +- ; 

The set (jlOp for ?/ > with Vn' at = is solved by matching in the interme- 
diate region the corresponding solutions ''^n^} one of which is valid in the 
inner region close to the impurity centre and the other represents a solution of 
the outer region distant from the centre. This method was originally developed 
by Hasegava and Howard [26^ in studies of excitons subject to strong magnetic 
fields and then successfully employed for the investigation of the impurity and 
exciton states in quantum wells [2 7) , super-lattices |28j and quantum wires |18j . 



3.1 Inner region 

In the inner region 

< y ^ ro {tq — p\e% — _£;2|-i/2 jg fj^^ effective size of the Coulomb state) an 
iteration procedure is performed. The first integration of the set (jlOp . in which 
we neglect the terms consisting of the energies e^r and E, with the trial functions 

(1) (1) (2) (2) . 

'"no = o-n y'"m = °-N gives 

4liy)^^'-N^~4^^idiFih)+d,Fih)] ; /i,2 = |^ (15) 

where 



/ \ 21/1 + 1) , |/|»1 

and where 

is the dimensionless strength of the Coulomb potential. 
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The function w^j (y) can be obtained from eq. (fT5|) by replacing a^^^ by a j^'' 
and a^'* by — a^'. Subsequent integration leads to the two independent partic- 
ular solutions {wjv+'''^w+} ^^'^ corresponding to the relationships 
a^'' = ibia^''. The linear combination of these solutions provides the general 
iteration functions, which read in the region y ^ di 2 



^Sn(y) = MQ+Cn) ; v'-^ljy) = Rn cos(Q+Cjv) ; Q{y) - g 1^ (in ^ + 1 

^ (16) 

where 

I? = yd;d^exp|^(di-d2)ln^| (17) 

and where Rn and are the arbitrary magnitude and phase, respectively. 
Since the potentials ([TT|) satisfy VjyN'iy) — VNN'{—y) the wave two- vectors 
vjv{w^\ w^''} are classified with respect to parity. Further we focus on the even 
wave vectors for which IIujv = vn where 11 = ttct^ with TTv{y) = v{—y) {gz is 
the Pauli matrix). The condition of even parity imposed on the wave vector 
^jvin formed by the components (jl6p . implies that the phases Ctv are equal to 
the half integer of tt. 

Obviously, an alternative way to derive eq. (fT6|) is to solve eqs. (fTO|) for 
V/VAr(y) ~ — £n = E = Q and then to compare the resulting solutions 

i'/^+^_(y) ~ Ay^^'^, expanded in series up to the terms of the first order of q, 
with those given by eqs. ([15]). The calculated constant A leads to the functions 

m- 



3.2 Outer region 

a) Discrete states 

The exact solutions to eqs. (fTO|) at Vn'N = ioi N N' in the region 
y 2> (ii,2 with VNN{y) — —Py~^ are calculated by the same method employed 
in studies of a relativistic electron in hydrogen [2 9) and super-heavy atoms with 
the nuclear charge number number Z > 137 [iTl [121 HSl [H] 



A 



N 



cosh '■ 
sinh%T-^ 



(18) 



where 
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T = i/Ary, z/AT = -a/^I^ - £:2; tanh?/) = K = 7? - i; ^ = iq; 



P^N 2T{r]N)^ cosh ^/^^r h + ^ tanh 

770 <C 1 is the quantum number labehng the ground state 

and where Wk.i_i{t) is the Whittaker function associated with the Kummer func- 
tion U EHl 



W^At) = e * r5+^;7(a, c, r); a = i + /i - k, c = 1 + 2/i, (19) 



with [21] 



T""*; r > 1 

(7(a,c,r)-<^ ^11^ , £(£^^i-c.Rec-l- r^l- r«l ^^^^ 
[ r(a-c+i) + r(a) ^ ,nec— i, c^i, r <^ i. 

The functions (|18p are normahzed to + W/^''^)'^?/ = 1- 

The asymptotic behavior of the outer functions (|18l) at large distances y ^ 
follows from eqs. (fT8)) . (fT9|) and (|20)) 

7;(,i'^)(y) ^ exp ( -]L + ^^^i^^Ji) ; ro = — . (21) 



In the region r <C 1 eqs. p8)) . (|19p and ([20| lead to the expressions 

^wout = '^^'^ '^^ ' "wout " ; WAr(T) q In r + 6Ar, (22) 

where 



On = argr(-2ig) +argr(-7?jv +i(?) -arctan^ ^^1 + ^ - ij , (23) 

and Pjv are arbitrary constants. 

Matching the eqs. ^TE\\ and in the overlapping intermediate region 
d <C y <C To we impose the condition 



Vin(2^) _ Vout(y) 

which yields 



""nIiM "Sout(2^) 



LON -Q-Cn = sTr; s = 0,±l,±2,.... (24) 
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Using the properties of the arguments of the F-functions in eq. ([^ for 
a small parameter g <C 1 and for the quantum numbers 77 at = n + SNn,n = 
0,1,2,.. .,(5 jvn < 1 [30] (see also Ref. [ITJ[T3] for details) and choosing — ■7t/2 
we arrive at the equation for the corrections Sj^n 



q 



arctan ■ 



arctan 



2{n + SNn) 



ln(n + 5n7i) 



+ V(l + n)+ In j'""^ + 2C - 1 ^ 0, (25) 



where n = 0, 1, 2, . . . , ^(1 + n) is the psi-function (logarithmic derivative of the 
F-function), C — 0.577 is the Euler constant. The corrections Snu calculated 
from eq. p5|) determine the Rydberg series of the discrete energy levels i?Ar„ 
adjacent to the size-quantized energy level En 



Enu = 



--N 



2(n+<5jv„)^ 



n= 1,2,.. 
n = 



(26) 



1+1 



which allow to estimate the size of the Coulomb state in eq. (pij) rg ~ d/\N 
a\'Kq for n = 1, 2, . . . and 



''0 



for n = 0. 



1+1 



Clearly from (|26|) . the existence of the intermediate matching region di^2 ^ 
J/ ^ ro is provided for excited state n = 1,2,... by the employed above small 
parameter q and for the ground state n = by the condition 

(1 ± ^)|iV - a|7rzo < 1, i.e. zq = < 1- 

a ojvo 

The correction 6n{) = q/ zq satisfies the transcendental equation 



In zq + - 



arctan zn — arctan — 
2 



m^^^ + c-i.o, ,27) 



while the corrections Snu for the excited states n ~ 1,2,... can be calculated 
in an explicit form 



In g + In n + — - V(l + ri) - In "^^""^ - 2C + 1 

2n 2d 

b) Continuous states 

Since our approach to determine the wave function of the continuous states 
closely resembles that applied above for the wave functions of the discrete states 
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only the basic points will be given below. Setting in eqs. ([T5| = — 2ifc we 
obtain 



it) 
it) 



B 



N 



cos^r 

sin 



w,At)~^'-^w,+,^,{t) 

W,,,it)+i'-^W,+,.,,{t) 



where 



1 / pk 



pk ^ . q 
— ; K — 1 

En sin ip 



q 



2tt tan ip 



exp 



irq 
simp 



(29) 



(30) 



The wave vectors ([M]) are normalized to S{k — k'). 

At large distances ky ^ \ the wave functions (j^ . , (^(7)) have the 

asymptotic form of the outgoing waves 



^%+^ iv) ^ exp <! ik^y + i^^ In 2kNy 



. q 
\— — 

sin ip 

Further we introduce the real wave functions associated with the standing 
waves 



Vout'^'^'' ~ ^ ^^jv+lij + e (t) 



where the functions v^j^'^\t) — v^j^'^^*{t) have the asymptotic form of the ingoing 
waves ~ exp |— ifcy — i^j^ In 2fcy| and Dn and are the arbitrary magnitude 
and phase, respectively. 

In the region \t\ < 1 eqs. ([291), (HH) and ^ lead to 



,(1) 

AfOut^ 



eli ^5^M cos (Hat + g In |t I + ^4 



e""^^ ^^Srn^ cos (-f^AT + gin |i| + e-) 



(31) 



where 



r+.- = r 



-1? 



± 1 



sm ip 



= argr(-2i(7)Targr+,_;M(^) = 



1 + sin (f 



(2,) _____ 

The wave functions ''^JyQ^^lO can be obtained from eq. (PT|) by replacing Z?jv 
by -Dn and ^+._ by ^+^_ - 7r/2. 

For q <^l, (p the wave functions (l5Tt read 



„(i) 



AfOUt i^^ ^ [^^"^ ~ "^"^ ' Vout (^^^ ^ t*^*^^ + '^^^ 

(32) 
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with 



XN = qin \t\ + + AN = nN- ^(e- - cn = 9^(1 + coth . 

^(33) 

Similar to the case of the discrete states we obtain the equation for the phase 
An on equating the ratios (y) / v'"^\y) taken for the iteration (fT6|) and outer 
functions ([5^ 

cot Aat = — tan(x - Q - Ca?)- (34) 

Cat 

Since [30] 

h- - e+) = ^ - f^^"^ ^y(-- arctan^) 

eq. (j34l) acquires for Cjv = 7r/2 an exphcit form 



cot Af = 









hcothf) 




In^ - i 




+ i^V 




-2C+1 



(35) 



Since at g <C 1 the phase fi^^^ 
with C(3) 

in eq. 17(0' = 0. 



|C(3)^ (el's) is the Riemann zeta function 



1.20) is the value of the higher order of smallness ^ q^ <^ 1 we set 



(1 2) 

As expected, setting in the functions v 



jvout 



kN (r+ = T*_) and 



(1 2) 

then matching these functions with the iteration functions v •', 
° wit 



16]) we obtain 

the equation for the phase ■ Substituting this result into equation cot rijv = i 
determining the poles of the S matrix with S = exp(2inAr) [TnUSH] we arrive at 
eqs. for the discrete energy levels. 



4 Double-subband approximation 

Below we consider the coupling between the ground iV = and first excited N = 
1 states described by the system of the four upper eqs. ([T^ at K_io — V-\\ = 0. 
Applying the iteration method with the trial functions Vq — and v\ = 

(12). 

a\ ' we arrive at two particular linear independent four-vectors, having the 
components v^^f^ calculated for a^^ ^ = iag^\ a[^^ = ia[^^ and = —ia^p, a^-* — 

—ia^i K The linear combination of these vectors taken for a^^ l — i?o,i exp [i(Co,i ~ 
provides the general expression for the iteration four- vector with the components 



"Oit 
„(2) 



(y) = Rq sm{Q + Co) + Riqioi cos Ci 
(y) = i?o cos((3 + Co) - Riqim sinCi 



(36) 
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where i?o,i and Co,i are an arbitrary constant and phase, respectively. The 
parameter 



7^70,1 = cos aa 
+ sin ao 



Ci(| + <.o)-Ci(|-ao) 



(f 



ao = 



(37) 



consisting of the integral sine Si and cosine Ci [30], describes the coupling in- 
duced by the potentials Vbi — Vio (ITlT) . The functions ^^|^(y) and w|j^(y) can 
be obtained from the functions v''^^ (y) and w^?"* (y) p6p by mutual replacing 

Equating the ratios of the functions of the continuous spectrum 
'i'oout(*)/"oout(^) ^^'^ iteration functions dSS]) ^^^^Jjj) / 'v'i^^r.iv) and then 



om^ 



matching the ratios of the functions of the discrete states 



„(i) 



,M/-inU.(-) 



and the iteration functions ''^{{■^{v) I ''^{{■^{v) obtain the set of equations 



Ro [sin(xo - Q - Co) - Co cot Aq cos(xo - Q - Co)] + 1 

Riqioi [cos(xo - Ci) - Co cot Ao sin(xo - Ci)] = ; > (38) 
Roqioi cos(iLJi - Co) - ^1 sin(i:ji - Q - Ci) = , J 

where the functions xn, Qt^n are defined by eqs. (155)) . p^ . (f^ . respectively. 
In the limiting case of negligible coupling the set (155)) decomposes into two 
independent equations relevant to the discrete (1^ . ([55]) and continuous ([Ml) 
states. 

Solving the set ()38p by the determinantal method we obtain the equation for 
cot Ao, which is then expanded in series with respect to the parameter q <^1. 
Keeping at Cat = 7r/2 the terms of the first order ^ q we arrive at the equation 
for the phase f2o in an explicit form 

^-o co(wi-g-f) 
'^"^ "° ^ (-r-Q-i)(0-Xo + f)-.^7o^/ 
Substituting cq. ^ into equation [TUII^IS^ 

cot rJo = i (40) 

the complex quantum numbers rji introduced in eq. ()18p can be calculated, 
which in turn determine the complex energy levels adjacent to the size- 
quantized first excited energy level ex — {1 — (j)TTp / d 

Em =£i-£i ^, f .2 i^; n = 0,l,2..., (41) 

2(n + oiJ 2 

where the second term in the right-hand part is the Rydberg series of the energy 
levels associated with the quasi-Coulomb diagonal potential Vii{y) ([H)) (no 
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coupling). The following notation in eq. (|1T|) for the resonant shift Win and 
resonant width Tin both induced by the inter-subband A'^ = 0, 1 interaction is 
used 

Tin = 2£l , ,3 G«('?ln)g70 1^0.1 (42) 

and 

^1" = J^! G„(^i„)g'7o'.i^o,i- (43) 



In eqs. (1121) and (gS]) 



Go-\'5io) = '5ro' + ('^?o+'7')"'-(2'5?o+29')"': G„(5i„) = <52„+g2, n = 1, 2, . . . , 

where the corrections Sin can be calculated from eqs. ([25]) ([ST)) ([28| . In the 
logarithmic approximation glng^^ <C 1, Go((5io) = 2(5Jo, G„((5i„) = (5^„, n = 
1,2,.... The quantum number feg can be found from equation 

with d?]) for £0,1. 

In conclusion of this section note that the equation absolutely identical to 

_____ fl 2) 

eqs. (PT|) - P5)) can be derived by matching the real iteration functions v^--^^ iu) 

fl 2l _____ 

and complex functions of the continuous states Uq^ (t) ([^^ having the 
asymptotic form of the outgoing wave. 



{n + Sin^ 



5 Three-subband approximation 

In this section we consider the coupling between the discrete states adjacent to 
the highest size-quantized level £_i and and the continuous states attributed 
to the low- lying levels Eq and ei. Below we neglect in the set the off- 
diagonal potentials Vbi and Viq describing the interactions of the = 0, 1 
subbands. Extending the iteration procedure employed above for the single- 
and double- subband approximation to the present stage with the trial func- 
tions Wjy' = aNi.2 we arrive at two particular linear independent sixfold vec- 
tors V-^^^{vq^\v^\v[^\v^^\v'"^1,v^^I) calculated for a^l,, _ — ±ia^l^ - , N = 
0,1,-1. Taking a^^\ _ = 

i?Arexp [±i(C7v — ■§)], we obtain the components w^'^'' of the total iteration 
sixfold vector ^ V+ + V- 
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v^^liy) = ^0 sin(Q + Co) + ^-1970-1 cosC-i ; 

vlfjy) = ^1 sin(Q + Ci) + RiQli.-i cosC_i ; } (44) 

lin^y) = ^-1 sin((3 + C-i) + ^0970,-1 cos Co + Riqji,-i cosCi , 

where Q{y) is determined in eq. (|16p and Rn and Cw are arbitrary constants 
and phases, respectively. The parameter 70,-1 — 70,1 (|37p . while 



27r7i,_i = cos2q:o [— Si(7r + 2q!o) — Si(7r — 2ao)] 
+ sin2ao [Ci(7r + 2ao) - Ci(7r - 2ao)] , "o = ^ (45) 

describes the coupling induced by the potentials V-n — V\-\ (fTTj) . The func- 
tions ^^1^-^(2/) can be obtained from the functions ^^^•|-(2/) (HUj respectively, by 
replacing sin((5 + Ca^) by cos((3 + Cw): cos(CAr) by sin(CAr) and g by — q. 

As mentioned above further we match the wave functions of the continuous 
spectrum v-^' iv), N = 0,1 having the asymptotic form of the outgoing wave 
Vn'+ (|29p to give in turn the sixfold vector V^^t with the components 

^wout (y) - (1 + Civ) cxp [i {q In 2A:jv2/ + M] 
^{l-CN)exp[-i{q\n2kNy + ^-)], N = 0,1, (46) 

where cn and a-re given in eqs. ([5^ pip and ([2^ . The wave functions 
^Nout^y"^' = 0, 1 can be obtained from the functions '^^^out^^-' (HUli respec- 

f 1 2) 

tively, by replacing by £,+.- + \ . The wave functions "^J-YQ-yy^hi) have the 



form (|22|) . in which j/at (llSp and Ojv psp are calculated for iV = — 1. 

Matching the sixfold wave vectors V\-^ and Vq^xX within the intermediate 
region by imposing the conditions 

jvout _ jvm AT — n 1 —1 

(2) (2) ' ' ' 

Noui Nm 

fi 2^ I t I t (1 2^ 

where wJ^q^^ are given by eqs. (|46| and ([TS]) and ^^J^-j^ by eqs. (f44|) we obtain 

i?o[4"^ cos(Q + Co) - i4^' sin(0 + Co)] - i?-ig7o,-i [4"' sinC-i + i^^^ cosC-i] = ; 
i?i[4"^cos(Q + Ci) -i4^'sin(Q + Ci)] -i?-ig7i,-i[4"^sinC-i +i4^^cosC-i] - ; } (47) 
sin(a;_i - Q - C-i) - Riqji-i cos(tJ_i - Q - Ci) - RmiQ-i cos(a;_i - Q - Co) = 



In eqs. (|47] 



^iv^' ^ = (1 + cn) exp [i (g In 2kNy + C+)] ± (1 " cw) exp [-i {q In 2fcAry + C-)] 
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and Q{y) and w_i(y) are introduced by eqs. ([TC]) and ([^ for = — 1, respec- 
tively. 

Solving the set (|T7|) by the determinantal method we obtain the equations 
for the complex quantum numbers rj^ = 2qE /pv^i (|18p . in which we take for 
the phases CAf = f, iV = 0, 1,— 1 and keep the terms of the first order oi q<^\ 

UJ-l-Q-^ = q^ 7Ar,_l(gr^jV,-l+iBN,-l), (48) 

N=0,1 

with 

- Bk-i (i- ^ + 1° 7 + ^ - l) ' = 

and the quantum numbers fco.i are obtained from 
=el+ p^kl =el+ p^kl = el 

with d?]) for e_i. 

The complex quantum numbers rj calculated from eq. (j48p lead to the com- 
plex impurity energy levels adjacent to the size-quantized second excited sub- 
band £-1 

i?_l„=£_l-£_l— -+W-^n-l^^. 71 = 0,1,2,... (49) 

where the resonant width r_i„ and shift have the form 

r-m - 2£_i f G„g(7iVi-gi.-i+7oVi-5o.-i) (50) 

[Tl -\- — In ) 

and 

W.^n = -e-i . JI ^^ Gnq' (7?_iAi,_i -)-7o-i^o,-i) ■ (51) 
The coefficients G„ are defined in eqs. and (|43p . 

6 Discussion 

We define the binding energy of the electron E^^^ in the n-th quasi-Coulomb 
state associated with the N size-quantized subband as the real part of the 
difference between the size-quantized energy e^r ([7]) of the free electron and the 
energy of the impurity electron E^n given by eqs. (pS)) . (HJ) and (|1T|) for the 
ground iV = 0, first TV = 1 and second TV = — 1 excited subbands, respectively. 
Since the resonant shifts Wf^n l(43)) . ((5T|) are of the order of ^ 1 (Wo„ = 0) 
with respect to the Rydberg energies determined by the second terms in the 
right-hand parts of eqs. (|¥T|) and (PT]) the binding energies read 



(n + 5_i„)2 
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d ^ (nm ^) 



Figure 2: The binding energy EqJ (|52p of the ground state n — calculated 
for((7 = 0.13; 0.24) as a function of the reciprocal width ^ of the GNR. Impurity 
is placed symmetrically to the boundaries (a;o = 0)- The parameter a = 0.3. 



Enu = 



En 



n = l,2,... 

n = 



(52) 



where the corrections Snh can be calculated in the single-subband approxima- 
tion from eqs. ^7} and ([^5]) for the ground n = and excited n — 1,2, .. . 
impurity states, respectively. It follows from eqs. (|52]) and ([7]) that the binding 
energy ^ en ^ and the osciUatory part of Sn (see a(d) in eq. ([7])) 
decrease with increasing the ribbon width d. In an effort to render our calcu- 
lations close to an experimental setup, we take below for the estimates of the 
expected values the parameters q — 0.13 (e ~ 25) [33] and q = 0.24 (e ~ 10) 
corresponding to the IIf02 and sapphire, respectively, employed as substrates 
for GNR [34^. The latter parameter q is close to the limit caused by the condi- 
tion Zq (see below ea. ip^ ). Further we focus on the monotonic dependence 
~ d~^ and keep the the parameter a in eq. ([7]) for the levels to be a ~ 0.3. 

The dependencies of the binding energies on the width of the GNR d for the 
ground state for different strengths of the impurity potential are given in Fig. 
2. These graphs, while ignoring the oscillations, are qualitatively completely in 
line with the data of the numerical calculations and experimental observations 
recently performed with the related Coulomb systems. The exciton effects in 
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4 0.8 



Figure 3: The dimensionless binding energy Ej^^ = Ej^yirpdr^ calculated from 
(O, (HH), ^ with q = 0.20 for the ground state {n = 0) plotted as a function 
of the effective index of the corresponding subband v = \N — d-\ and the relative 
impurity position s = ^ in the GNR of width d. 



the armchair GNRs were studied in frame of the tight-binding model [35j and 
density functional theory [36], while Han et.al. [34] investigated experimentally 
the influence of the localized states in GNRs on the electron transport. The 
relation E'-'^^ ^ d^^ including oscillations [35] [36] have been found to occur. 
The differences between the impurity states considered here and the exciton 
and localized states prevent us from a detailed quantitative comparison. 

The Coulomb pattern of the energy levels ([26]) enables to introduce the 
effective Rydberg constant Ri/n, the Bohr radius aoN and the mass AIn for the 
impurity electron in GNR 

q^lN-alirp d _\N^a\7^ 

2d ■k\N — a\q pd 

which additionally illustrate the physical reason of the bonding of the impu- 
rity electron, namely the quasi-lD geometry of the GNR. Note that the bound 
states arise at any finite width d < oo. This result is qualitatively analo- 
gous to the effect of anti-diamagnetism caused by the influence of the magnetic 
field on the weakly bound atomic state. Demkov and Drukarev [19, consid- 
ered the 3D potential well of small radius and depth to provide the capturing 
of the electron. It was shown that the arbitrarily weak magnetic field B in- 
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duces the bound electron state with the binding energy 
ifi/eBY^'^ is the magnetic length). The common reason for this is that the 
finite width d < oo and magnetic length as < oo transform the graphene mono- 
layer and atomic structure, respectively, into the quasi-lD systems, which are 
more favorable to generate bound states. The dependencies E'^'^^ ^ and 
£;(b) correspond to the different dispersion laws namely i?'-^' ~ and 

£;('') p (jp ^ fijT is the momentum) for the atomic (r ~ as) and GNR [r ~ d) 
electron, respectively. 

The dependence of the binding energy (|5^ on the displacement of 

the impurity centre xq from the mid-point of the ribbon a; = is contained 
in the corrections (5Ar„(xo) namely in the term InDd"^ in eqs. ((27|) and ((28l) . 
while e'^^^ as a function of the effective number of the subband |iV — crj is 
given by the sub-band threshold ejss ^\ mainly and the term In | — (t| in the 
correction (5Ar„. The dimensionless binding energy E^Q/Trpd~^ as a function of 
the effective quantum number |iV — a\ and relative displacement xo/{d/2) for 
the ground n = state is depicted in Fig. 3. Clearly, the higher the subband 
i.e. the greater the value |7V — a\ is the less the binding energy i^^g. Also 
the binding energy decreases when the impurity shifts from the mid-point of 
the ribbon towards the boundaries. The latter conclusion coincides with those 
obtained for the quantum well in Refs. [371 [Ml [Ml iOl EI] • 

The inter-band coupling shifts the strictly discrete excited Rydberg series 
Enji {N 7^ 0) ([26)1 calculated in single-subband approximation by an amount 
Wnu (dSj) N = 1, and ([5T|l = — 1 and transforms them to the quasi-discrete 
levels of width Tnti ([12|) = 1, and ([Mj) iV = -1. Note that the conclusions 
made on the base of the first and second excited subbands can be qualitatively 
extended to others. Since the resonant shifts Wnu ~ first are much less 
than the resonant widths rAr„ ~ gr at g <C 1 {Wnu ^ ^Nn) and second the 
resonant shifts do not change the discrete character of the energy spectrum p6p 
we focus on the widths Tnu- It is clear from eqs. ([42]) . and (|50p that the widths 
^Nn ^ £n ^ increase with decreasing the ribbons width d. Note that this 
dependence is opposite to that in a semiconductor narrow quantum well: the 
narrower the well is the less are the resonant widths [42l [38l [22] ■ The reason for 



this is that in the quantum well the resonant width Tj\in ~ ^nIi {-^Nn/^^N 



where the impurity Rydberg constant Ry ~ Ej^^ and the binding energy Ej^^ do 
not depend on the well width d, while the inte-rband energy distance Aen ~ d~'^ 
increases and consequently the resonant width decreases with the narrowing of 
the quantum well. For the ribbon E^^]^ ^ en ^ Aejv ''^ d~^ ([52|l and the inter- 
subband coupling do not depend on the ribbon width d and Tjyn ^ ^nIi ^ 



The dependence of the resonant widths ri„(a;o) (|42|). calculated in the 
double-subband approximation, on the position of the impurity centre xq is 
described by the coupling parameter 701 (137]) and the corrections Sin |27|^ and 
([^5)1 . The dependencies of the relative resonant widths Fin/ei ([1^ . on the di- 
mensionless shift s — 2xo/d for the ground n = state are presented in Fig. 4, 




En ^ d 
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Figure 4: The resonant width Fio (|42|l of the ground impurity state {n — 0) 
relatively to the corresponding threshold ei ^ (Fio = Fio/ei) versus the rela- 
tive impurity position s ~ in the GNR of width d providing the parameter 
a = 0.3. The parameter q is taken to be g = 0.13; 0.24. 



in which the limitation on the parameter s are caused by the condition imposed 
on zq placed below eq. (|26|) . For the impurity positioned at the mid-point of the 
ribbon xq = the resonant width and shift both vanish (Fi„(0) = W"i„(0) = 0) 
because of the even x-parity of the Coulomb potential V{p) ^ in eq. 
and opposite parities of the neighboring = 0, 1 transverse a;-states to give 
M31 = 7o,i = 0- Both in the quantum well and in the ribbon the shift of the 
impurities from their mid-points eliminates the even x-parity of the potential 
V{p} ^ in eq. that leads to the coupling between the N — and = 1 

subbands. If the impurity displaces from the mid-point towards the boundaries 
|a;o| = d/2 the resonant widths Fi„ monotonically increases. This corre- 
lates completely with the analogous dependence found for the impurity states 
in the semiconductor quantum well [42l |38l [27]. For small shifts ao ^ 1 in 
eq. ([37|l we obtain for the parameter 70,1 in eqs. ((42|) 70,1 = (f )'S'«(f )ao ^ 1 
while for the impurity positioned close to the ribbon edge a;o ~ | we obtain 
701 = ^Si{n) with Si{^) = 1.37, ^^(Tr) = 1.85. Note that the zeroth width of 
the first excited {N — 1) n— series in case of the symmetrical xq = impurity 
position is a consequence of the double-subband approximation. In the multi- 
subband approximation the levels of the above mentioned series would acquire 
finite widths. 

The dependence of the resonant widths (f50| of the impurity states cor- 
responding to the second excited subband = — 1 on the position of the 
impurity centre is completely different from that related to the first excited 
subband = 1. Equation ([50| shows that contributions to the resonant 
widths F_i„ arc caused by the coupling with the subbands A^ = 7o-i) 
([37]) and A^ = 1 7i,-i) (|45|). Note that the estimated contribution to 
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\s\ 

Figure 5: The resonant width r„io of the ground impurity state (n = 0) 
relatively to the corresponding threshold e_i (O (r_io = r_io/e_i) versus the 
relative impurity position 2|2. in the GNR of width d using the parameter value 
a = 0.3. The parameter q is taken to be g = 0.13; 0.24. 

the resonant width r_i„ caused by the neglected coupling between the N = 
and N = 1 subbands is of the order of qjQi ^ 1. In the vicinity of the 
mid-point (|a;o| <C d/2, ao <C 1) the subband = 1 contributes mostly 
(70, -1 ~ 0, 71,-1 — — i5i(7r)), while for |xo| ~ d/2, ao — 7r/2 both subbands 
contribute 70,-1 — ■^Si{Tr), 71,-1 ~ ^Si{2Tr). The position of the impurity 
xo — ^o^o, at which the effects of the subbands = and = 1 on the reso- 
nant width r„i„ are in balance is determined by the root ao of the equation 

if-iBi-i = 7o,_iSo-i. 

to give the result ao = 0.58, xq = 0.37|. The coupling between the subbands 
= — 1 and A^ = 1 provides the nonzero widths r_i„ and shifts W^in for any 
positions xo of the impurity. The width F-io as a function of the impurity shift 
Xo is given in Fig. 5 demonstrating the monotonic drop within the same regions 
as those corresponding to Fig. 4. 

As mentioned above the presented method is valid under the conditions 
q ^ 1 for the excited impurity states n — 1,2,... and ^0(9) ^ 1 (|27|) for 
the ground state n = 0. Under these conditions the radius of the impurity 
state considerably exceeds the width of the GNR so that the ribbon is narrow 
compared to the impurity size. It follows from eq. (|25|) that a "big logarithm" 
can only hardly achieved [TTJ [T2j [131 E] i-e. the logarithmic approximation 
jlngl ^ 1 ensures the real smallness of zq. However the previous calculations 
related to the ground state of the quasi-lD diamagnetic exciton [13] and present 
estimates show that a reasonably small parameter q leads to values zo < which 
provide a quite accurate and adequate description of the ground impurity state 
in GNR. 

Taking into account possible experiments we estimate the expected electron 
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Figure 6: The dimensionless binding energy EqJ = Eqq /eq of the quasi-classical 
ground state {N = n = 0) found from ([55]) for Eqo and from ([7]) with a ~ 0.3 
for £o versus the parameter q. 



binding energy for the impurity centre placed at the middle point of the GNR 
of width 1 nm on the sapphire substrate as i^Qp"* ~ 160 meV and on the Hf02 
substrate as Eq''q ~ 68 meV. This is less than the data attributed to the Si02 
substrate (e = 3.9) because of the relatively small screening of the impurity 
potential. Also an estimate of the lifetimes Tj^n — ^/^Nn yields for the impurity 
positioned at the mid-point of the GNR t_io ~ 0.21 ps and 0.049 ps for the 
Hf02 and sapphire substrate, respectively. For the Si02 substrate the screening 
of the impurity attraction is less, the lifetime is reduced and therefore less 
favourable for a corresponding experimental observation. A shift of the impurity 
centre \xo\ — OAd/2 generates lifetimes tiq of the same order as r_io at xq = 
0. The electrons captured onto such short-lived trap states will most likely 
contribute to the dc transport. However, the high-frequency response of such 
electrons may reveal the signatures of localization. 

Clearly, the above considered quasi- Rydberg series (|26|) do not cover the total 
set of discrete states. The oscillations of the wave functions and ([^ caused 
by the logarithmic term are an indicator of additional energy levels positioned 
below the series ((26|) . Since the possible strong shift of these levels away from 
the threshold £jv is against the spirit of the employed adiabatic perturbation 
theory {q <^ 1) describing the shallow energy levels we are limited to qualitative 
estimates based on the quasi-classical relativistic approach [9^ and [44^. 

In an effort to elucidate the origin of this additional series let us consider 
the so called "logarithmic" energy levels governed by the logarithmic potential 
taken for di = d2 = d (xq = 0). These levels can be calculated from the 
Bohr-Sommerfeld quantization rule 

rvo I 

riy)dy^Trh{n+-); 0,1,2,..., (53) 
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where V{y) and yo d are the classical relativistic momentum and turning 
point, respectively, with 



1 

7^ 



-N 



V\y) : 

Vp \ a a J 

Viyo) = 0. (54) 
Equation ([55)1 admits an exact solution which provides for the energies 



E 



Nr. 



2p 



In 



in 



\N- 



In Ki{sn) 



sn 



ENd 

2^ 



\N- 



2q 



(55) 



where Ki{s) is the modified Bessel function [30 . 



The binding energy E 



Nn 



Nr. 



with E]sin calculated from (|55|) reads 
E^l^\ ~ q\a.q both for g < 1 (sa, > 1) and for g > 1 {sn < 1, ling] » 1). 
It follows that the weakness of the logarithmic singularity and smallness of the 
strength of the impurity potential (g ^ 1) seems not to provide the bonding 
of the quasi-classical relativistic electron {E^^^ < 0), while a sufficiently strong 

attraction {q < 1) could produce a localized impurity state (£^^)j > 0). The de- 
pendence of the binding energy of the quasi-classical ground state {N = n = Q) 
found from ()55p on the parameter q is depicted in Fig. 6. The ground "logarith- 
mic" level arises at the critical value qo ~ 0.48 and shifts towards lower energies 
to provide for the binding energy 0.1 < E'qqV^o < 0.5 for 0.54 < q < 0.74. 
The above can be considered as no more than only a qualitative evidence of 
existence of such additional states in GNR that have transformed from the col- 
lapsed states in the graphene monolayer governed by the 2D impurity potential 
~ —r~^ Though the "logarithmic" and quasi- Rydberg levels in principle 
correspond to the same region of the parameter q < 1 any numerical comparison 
between the results for the quasi-Rydberg series based on the Dirac equation 
and those for the "logarithmic" levels derived from the quasi-classical method 
applied moreover to the ground state seems to be incorrect. The total set of the 
impurity states in GNR requires a further study of the equations ([TU)) with the 
potential (llip , having the logarithmic singularity in the vicinity of the impurity 
centre. 



7 Conclusion 

We have developed an analytical adiabatic approach to the problem of bound 
and meta-stable (Fano resonances) quasi-Coulomb impurity states in a nar- 
row gaped armchair graphene nanoribbon (GNR). The width of the GNR is 
taken to be much less than the radius of the impurity. This adiabatic criterion 
implies a variable width of the GNR and simultaneously the smallness of the 
Coulomb interaction relative to the size-quantized energy induced by the GNR. 
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The energy spectrum of the impurity electron is a sequence of the scries of the 
quasi- Rydberg discrete and resonant states adjacent to the ground and excited 
size-quantized subbands, respectively. The binding energies and the resonant 
widths and shifts attributed to the inter-subband coupling are calculated in 
an explicit form in the single- and multi-subband approximation, respectively. 
The binding energies and the resonant widths both increase with decreasing the 
GNR width. As the impurity centre displaces from the mid-point of the GNR 
the binding energies decrease, while the resonant widths of the quasi-Rydberg 
series associated with the first/second excited sub-bands increase/decrease, re- 
spectively. Our analytical results are in complete agreement with those found 
by other theoretical approaches and in particular numerical studies. Estimates 
of the expected values show that the bound and meta-stable impurity states in 
GNR can be observed experimentally. 
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